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Abstract. In this paper, we provide a probabilistic interpretation of the Volkenborn integral; 
this allows us to extend results by T. Kim et al about sums of Euler numbers to sums of Bernoulli 
numbers. We also obtain a probabilistic representation of the multidimensional Volkenborn 
integral which allows us to derive a multivariate version of Raabe's multiplication theorem for 
the higher-order Bernoulli and Euler polynomials. 



1. Introduction 



The Volkenborn integral was introduced in 1971 by A. Volkenborn in his PhD dissertation and 
subsequently in the set of twin papers [1]; a more recent treatment of the subject can be found in 
[S]. The Volkenborn integral, or fermionic p— adic q— integral on Z p , of a function / is defined as 



f 1+ % 

/ / (V) d^-q (y) = lim N V / (x) (-?) 

h v N^+oo l + q p ^ 



In particular, the q = 1-Volkenborn integral satisfies [2 eq. (1.6)] 



x) 



r 9 +°° 

where E n (x) is the Euler polynomial of degree n. 

Another interesting case is the q = 0- Volkenborn integral which satisfies [5, p. 271] 

/ e^d^{y) = f- i =Y. B ^)-< 



n=0 

where B n (x) is the Bernoulli polynomial of degree n. From this result we deduce the following 

Theorem. // / (x) is analytic in a neighborhood of then its q = 0- Volkenborn integral can be 
computed as the expectation^ 

(1.1) f f(x)dfi (x)=-Ef(x + iL B -l 



2 

where the random variable Lb follows the logistic distribution with density 
X ' (1-2) ^-sech 2 (ttx) , x € M. 

Moreover, its q = 1- Volkenborn integral coincides with the expectation 

f (x) d^ii (x) = Ef (x + iL E - i 

where the random variable Le follows the hyperbolic secant distribution with density 
(1.3) sech(irx) ,iel. 

Proof. These results can be proved by computing the characteristic functions associated to the 
logistic distribution 

Eexp (itLs) = -csch 



^in the rest of this paper, we use the notation Ex/ (X) for the probabilistic expectation f f (x)px [x) dx where 
px is the probability density function of the random variable X. When no ambiguity occurs, we also denote 
E/ (Xi , A~2 , • • ■ ) the expectation over all random variables that appear as arguments of the function /. 

1 
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and to the hyperbolic secant distribution 

Kexp (itL E ) = sech ^— 

□ 

The special case / (x) = x n yields the following moment representation for the Bernoulli poly- 
nomial of degree n 

(1.4) B n (x) = E (x + iL B - i 

and 

1 



E n (x) = E \x + iL E - 2 

for the Euler polynomial of degree n. Moreover, choosing x — yields the following moment 
representation for the n— th Bernoulli number 

B n = B n (0) = E (iL B - i 

and 

(1.5) E n = E n (0) = E (%L E - X - 

for the n— th Euler number where Lb and L E follow respectively the logistic distribution (11. 2p 
and the hyperbolic secant distribution (|1.3[) . 

An important feature of the logistic and hyperbolic secant random variables is the following 
cancellation property: 

Lemma 1. IfU E is a continuous random variable uniformly distributed over [0, 1] and independent 
of Lb then 

E (x + iL B - i + Ub) =x n . 

Accordingly, ifU E is a Rademacher distributed random variable (Pi {U E = 0} = Pr {U E = 1} = ^) 
then 

E (x + iL E - i + U E 

Both results can be shown considering the characteristic functions of the involved variables; for 
example, in the Bernoulli case, 

Eexp (t (iL B -^\\ exp (tU B ) = 1 

so that all integer nonzero moments of the random variable iL E — \ + U E vanish. 

The Volkenborn integrals were used by Kim et al [lj to obtain non-trivial identities on Euler 
numbers E n using integrals of the Bernstein polynomials defined as 



B k ,n (a) = y k )x* (1 - z) , < x < 1. 

In this paper, we show that the probabilistic representation (jl.ip of the Volkenborn integral 
makes its computation very easy to handle. We illustrate this fact by extending the non-trivial 
identities of pQ to the case of Bernoulli numbers. In the second section, we derive the probabilistic 
equivalent of the multidimensional Volkenborn integrals as introduced in [5] and we use it to prove 
a multivariate version of Raabe's multiplication theorem for Bernoulli and Euler polynomials. 



2 note that it differs from the n— th Euler number of the first kind defined by E n = 2 n E n (i) 
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2. Identities for Bernoulli numbers and polynomials 

2.1. first-order identity. In order to obtain non-trivial identities on Bernoulli numbers, we re- 
place the Bernstein polynomials used by Kim et al by the Beta polynomials 

*8 fe ,„ (x) = x k (1 + x) n ~ k , < k < n 

and, with X = iLb — \ 1 compute the expectation EQ3fc >n (X) in two different ways: 
- the first way is by applying the binomial formula 



'n — k\ • r i—^ (n — k 



(2.1) E<B k , n (X)=EX k J2 ■ P = E ' ) B 



j=0 " 3 ' j=0 



1j+k 



- the second way is by expressing X = (X + 1) — 1 so that 
(2.2) E<8„, fe (X) = e£ ( k ) (1 + Xt- k ^ k ^ = £ (*) (-1)' {S„_, + 5 n _,_ x } ; 

since j < k < n, the Kronecker adds a term ("l}) (—1)™ 1 if fc = n — I or ( "i) ( — 1)" 1 if = 
and no term otherwise. We conclude the following 

Theorem 2. The Bernoulli numbers satisfy 



Elo (5) M)' R n-j if0<k<n-2, n>2 



J2 ( n k ) B j+k = \ 2 = o (f) Bi] + (-1)- 1 k = n - 1, n > 1 



We remark that the case k = n reads 

n—l 



S " = EL (-l)^n-,+n(-l) 

2.2. polynomial identities. From p.4p . we deduce that the above results can be extended to 
the case of Bernoulli polynomials by choosing X = x + iLb — h ■ We deduce the following 

Theorem 3. The Bernoulli polynomials satisfy, for all < k < n, 

(2.3) J2 ( n . k ) B j+k (x) = E ( k ) (-I) 3 Bn-i (a) + - (n - fe)) x^- 1 (a - l) fc_1 

Proof. The left-hand side is a direct consequence of that of (|2.ip : the left-hand side of (|2.2[) with 
X = x + iLb — \ yields 

W*f) (-1)' 5 n -i(* + l) 



and since -B n -j (x + i7) = x™ J ', we deduce B n _j (x + 1) — £> n -j (x) = (ft — j) i n J ' 1 . replacing 
in the above sum yields the result. □ 



We note that the case k = reads 

'n 



Y J [ ] )B 3 {x)=B n { 



x) + nx n \ 



which can be restated as 

B n (x + 1)- B n (x) = nx"- 1 
and is nothing but the expression of the cancellation principle 

ES„_ X (x + U B ) = x n -\ 
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3. A POLYNOMIAL EXTENSION TO KlM'S IDENTITY 

In PP, the following identity is derived using the Bernstein polynomials Bk,n i x ) — (fe) xfe (1 — x) n ~ k 
E (" T k ) (-D J E j+k = ± ('-) {-!)*-> E n . A + 24. 

We now provide the following polynomial extension of this identity 
Theorem 4. T7ie Euler polynomials satisfy 

E ( n : fc ) (-i) j ^ (*) = (-i) n+fc+1 E («) + 2xk a - *r fe • 

j=0 \ •> ' j = Q V?/ 

Proof. We start from the identity 



j=0 v J ' j=0 



obtained by expanding either (left-hand side) the (1 — x) n k term of the (right-hand side) x k 
(x - 1 + 1) in the expression of the Bernstein polynomial. Replacing the variable x by x 



X + %Le — \ and remarking that 

(1 - x) n ~ j =E [1 — X — [ iLe — 7.) ) = En-i (X - 1) 



2 

2 

with, by the cancellation principle, 

E n _j (X-l) + E n _i (X) = 2 (X - 1) 

we deduce 



E ( ■) (-!) fe ^ (! - = E (- 1 )' £ ^ {~ E -i ( X ) + 2 - !)"^} 

= (-i) n+k+1 E (*) K-j (X) + 2J2 ( k ) (i - x) n ~ j 

this last sum being equal to 2X k (1 — X) n k , hence the result. □ 
We notice that the case X = is 

E ( n 7 k ) (-i)' = (-D" +fc+1 E ( ■) ^ + 

3=0 \ J / 3=0 VJ/ 



It can be shown that 



as follows: using the moment representation (jl.5l) the right-hand side reads 

E ^L E -l) n \l- (*L E - = E (_iL a - I)" * (l - (- t L E - ± 

by the symmetry of the hyperbolic secant distribution, and is thus equal to 

'L — k/ \ k 



with / (x) = (— 1)" i" fc (a; + 1) . By the cancellation principle 

Ef ( (iLe - 1) + l) + f (iL E - I) =2/(0) = 



so that the right-hand side if equal to — E/ (iLe — ^) which coincides with the left-hand side; 
hence we recover Kim's identity. 
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4. Multidimensional Volkenborn integral 
4.1. Introduction. In [2J, a multivariate version of the Volkenborn integral is denned as 

/ (x) d[i (x) = / ... / / (xi, . . . , x k ) dfin (xi) . . . duo (%k) ■ 



In particular, it satisfies, with y £ R fc and the notation |y| = X)i=i Vii 

This multivariate version of the Volkenborn integral can again be expressed as an expectation over 
a simple random variable as shown now. 

4.2. Moment representation and elementary properties. The Bernoulli polynomials (x|a) 
of order k and degree n with ifl with parameter a S M fc , also called Norlund polynomials, are 
defined by the generating function [3] 1.13.1] 

71 = j = l V 7 

and the corresponding Bernoulli numbers £?n (a) by 

(a) = (0|a) = fl 

In particular, taking a,j — 1 for all j £ [1, fc] and denoting 

B« (!) = flW (^|1,1,...,1) 

we deduce 



/■ + °° f- 

/ e ^l)^ (y) = E^ fc) W-r 



71=0 

We provide a multidimensional extension of the moment representation (|1.4p as follows 



Theorem 5. TTie Bernoulli polynomials B„ (Msl) satisfy 
(4.1) 



4*> (x|a) = E ^ + E a ^ - 5) j 



where the random variables I Lf) \ are independent and follow the logistic distribution (11.21) . 

1 J l<j<k 

As a consequence, the Bernoulli numbers Bn (a) satisfy 



(4.2) 



and £/ie multivariate Volkenborn integral, with x S K , 

& 



^ / (|X|) O> (X) = E / f E + lL B - \ 



Proof. Let us compute the generating function 



n=0 \ i=l V J V •' =1 / J =1 



with, for a logistic distributed random variable Lj, 

... (itU) ] 

El 

hence the result. □ 



A PROBABILISTIC INTERPRETATION OF THE VOLKENBORN INTEGRAL 



6 



The moment representation (14. ip allows to recover easily some well-known results about the 
higher-order Bernoulli polynomials. 

Proposition 6. The higher-order Bernoulli polynomials satisfy the identities 

(4.3) B« (x\a) = a n B n (£) , 

(4.4) £ ^ja'BS (i/|a) = B™ {x + „| a ) 
and 

(4-5) ^ +1 f ai + -; +flfc |aUo. 



Proof. Identity (|4.3p is a direct consequence of the moment representation (|4.ip ; identity (|4.4p is 
obtained using a binomial expansion of (|4.ip and identity f|4. 5[) by computing 

/ k \ 2n+1 



g 2n+l 9 " l a = E ' E a ' £ 



J B 

1=1 



and using the fact that the logistic density (ll.2p is an even function. □ 

We also deduce straightforwardly from a multinomial expansion of the representations (|4.ip and 
(j4~2| the following 



Proposition 7. TTie higher-order Bernoulli polynomials satisfy 

(x l + --- + x k \&)= E [■ n . )B ll {x 1 \a 1 )...B lk (x k \a k ) 
ii+...+i*=n \*i'-"'W 

and £/ie higher- order Bernoulli numbers 

Si fc) (a)- E (■ " • )B il (ai)...B ik (a k ) 

41 H Nfc — n 

These results extend Corollary 5 and Corollary 6 in [2] which correspond to the case a = 

4.3. Kim's identity for Norlund polynomials. In order to highlight the efficiency of the mo- 
ment representation (|4.ip . we derive now an extension of Kim's identity (|2.3p to the case of Norlund 
polynomials as follows. 

Theorem 8. For peN and < k < n, 

E ( n ; fc ) B $ k w = e f •) w + (» - i) • 

Proof. We start from the identity 



y& fc\ _ _ « _i_ t, x ^ i ^ 



n-j 



(4.6) e ; F j+fc =E ■ 

j=0 j=0 

and replace X by i + f^s^ — |J so that, from 14. ll the left-hand side reads 



i—k / j 
n — fe 



E( 7 K7w 



while the right-hand side is 



k 'A 



J'=0 
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Since ¥.B^_^_ 1 (x + U) = B^Zjly ( x ) with U uniform on [0,1] , we deduce by the cancellation 
principle 



77 -j 

which yields the final result. □ 

4.4. Kim's identity extended to multidimensional Euler polynomials. We now provide a 
multidimensional version of the polynomial Kim identity derived in Theorem [4] as follows: 

Theorem 9. The multidimensional Euler polynomials satisfy the identity 



3=0 x J ' j=0 

k 



+2(-ir +fe E(-)^- 

A— n \J / 



(x-l) 



Proof. Starting again from identity (|4.6[) , we take x = X + Ya=i {^^e ~ 5) and copute 

n—j 

»-3 p(p) 



(p \ 
l-x-E^-s)] 



= (_l)«-^^.(X-l) 



However, by the cancellation rule 

(X-l) + (X) = 22?£f (X - 1) 

and the result follows. □ 

4.5. Raabe's and Nielsen's multiplication theorem for Norlund polynomials. Raabe's 
usual multiplication theorem 

m — 1 



m 1 n B n (mi) = B n [ x 
1=0 



m 



and 



m_1 / I \ 
E n (mi) = (—1) E n (x H J , m odd 

/— n V 777/ 



i=0 

and Nielsen's multiplication theorem 

2 m_1 / / \ 

m~ n E n (mi) = 7 (— 1) B n+ i [ x H ) , m even 

n + 1 z — ' \ to / 

are an interesting feature of the Bernoulli polynomials since, as noted by Nielsen, [8, p. 54] 

It is very curious, it seems to me, that there exist polynomials, with arbitrary 
degree, that satisfy equations of the above form. However, it is easy to prove that, 
up to an arbitrary constant factor, the B n (x) and E n (x) are the only polynomials 
that satisfy the mentioned property. 
Using the moment representation and basic results from probability theory, we propose the follow- 
ing extension of Raabe's celebrated multiplication theorem to the multivariate case. 

Theorem 10. If men, 

m— 1 / -j k \ 

(4.7) m k - n B^ (mx|a) = £ B« U+-£aA|a 

h,...,l k =0 \ m i=l / 

and if moreover m is odd, 

m— 1 / \ 

(4.8) m-fjW (mx|a) = £ (-l) h+ - +h E™ a? +_£ 0< I < |a . 

ii,...,i fe =0 V m i=l / 
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Proof. Let us denote < U^' > a set of k discrete random variables independent and uniformly 

I J i<i<fc 

distributed in the set {0, ...,m — l\ and -J f/L \ a set of k continuous random variables 

independent and uniformly distributed on the interval [0,1] . For the Bernoulli case, we have 



toA * — ' \ to * — ' / \ * — ' V £ I to 

tl , — ,£fe=0 \ 2—1 / \ i—l i— 1 

to' 1 



71 \ i=l ^ ' i=l / 

= _L E ( mx + m £ ai - ±) + £ + E a, - 1) + £ rf) 

'' \ t=l ^ ' i=l 2=1 ^ ' i=l / 

Now we use the fact that 17 W + 17$ has the same distribution as mU$ so that we obtain 



m" 

\ 2=1 v 7 2=1 

and applying the cancellation principle, we deduce 

k 



For the Euler case, we need to use a signed measure (and then depart temporarily from the 

probabilistic context) defining the set \ > of k discrete variables independent such as 

I J i<i<fc 

each takes values in {0, . . . , k, . . . , to — 1} with a weig ht . Then 



k \ / k / ^ \ ^ k 



E (-i) h+ - +lk E n L + 1 £ =eu+^ ( ^ - 5 ) + - E 

ii,...,i fe =0 V i=l / \ i=l ^ ' 2=1 

2=1 ^ ' i=l / 



— E ( mi + m 



to" v 

= to" e f mx + m E to - 1) + E + E »* - 1) + E a ^ 

' \ i=l ^ ' i=l t=l ^ ' i=l 

where now are independent Rademacher random variables and, since to is odd, each U j 

has the same distribution as to U$ so that we obtain 

— E [ TOx + TO^a, (iL® ~ ^) + E° l \ 1 ^ ( e - 9 ) + m E fllE/ £ 



(0 



TO . . _ , _ , , 

\ 2=1 V ' 2=1 V 7 2=1 / 

and from the cancellation principle, we deduce the result. □ 

Raabe's identity (|4.7p and (|4.8|) are in fact given without proof in [71 eq. (1.6) and (1.7)]. The 
case for to even is not provided, so we prove now 

Theorem 11. With n £ N and to even, 

m—l / k 



~~k — r 
TO 



7 V ^' \2=1 / Jl,.».ifc=0 V i=1 ) 
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Proof. Let us define a variable W = {0, . . . , m — 1} with weights (—1)' . Then the right-hand side 
reads, with E>i = ilSg — ^ and Ei = iL^ — |, 




(k 
X + a * ( W i + mB * + E i + U e) 

5>i (Wi + mBt + ug) J 

i=l / 




but it can be checked that each Wi + takes values and m with respective weights \ and — | 
so that we obtain 



fe fc-i 



im" n (£„ ( x + mJ2a l B l + ^{Wi + Uf^j + ma k \ - E n (x + mJ^^B. 




— m 
2 

\ \ i=l i=l 

which coincides with 

^ / fc fe-i \ j / fc-i 

--m-"£7 n x + mj] a.B, + ^ (w^ + C^ ) + ma k U% } = - -m"" (ma fe ) x + ^ a t (w t + mB, + U 

\ i=i i=i / V i=i 

since makBk and ma,kU B cancel out. We are now back, up to a factor, to the same quantity as 
before except that n is replaced by n — 1 and fc by k — 1, hence the result. □ 
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